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Abstract. 

Recently, the QCD renormalization group (RG) equation at higher orders in MS-like 
renormalization schemes has been solved for the running coupling as a series expansion 
in powers of the exact 2-loop order coupling. In this work, we prove that the power 
series converges to all orders in perturbation theory. Solving the RG equation at 
higher orders, we determine the running coupling as an implicit function of the 2-loop 
order running coupling. Then we analyze the singularity structure of the higher order 
coupling in the complex 2-loop coupling plane. This enables us to calculate the radii 
of convergence of the series solutions at the 3- and 4-loop orders as a function of the 
number of quark flavours rif. In parallel, we discuss in some detail the singularity 
structure of the MS coupling at the 3- and 4-loops in the complex momentum squared 
plane for 0 < rif < 16. The correspondence between the singularity structure of the 
running coupling in the complex momentum squared plane and the convergence radius 
of the series solution is established. For sufficiently large n{ values, we find that the 
series converges for all values of the momentum squared variable = —q^ > 0. 

For lower values of rif, in the MS scheme, we determine the minimal value of the 
momentum squared above which the series converges. We study properties of 

the non-power series corresponding to the presented power series solution in the QCD 
Analytic Perturbation Theory approach of Shirkov and Solovtsov. The Euclidean and 
Minkowskian versions of the non-power series are found to be uniformly convergent 
over whole ranges of the corresponding momentum squared variables. 


PAGS numbers: ll.lO.Hi, 12.38.Aw, 12.38.Bx, 02.30.Mv 


1. Introduction 


It is known that the QCD rnnning coupling at the 2-loop order in MS-like (massless) 
renormalization schemes can be solved explicitly as a function of the scale in terms of 
the Lambert W function piaEi. The Lambert W function is the multivalued solution 
of 

Wk{z) exp{lTfc( 2 ;)} = 2 ;, (1) 

the branches of W are denoted Wk{z), k = 0, ±1,.... An exhaustive review of the 
Lambert-lV function may be found in ref. |1]. The relevant branch of W{z) which is 
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used to determine the coupling depends on the number of light quark flavours Uf. For a 
real positive momentum squared f (and above the Landau singularity if 0 < Uf < 8) 
the 2-loop coupling takes the form 




-(/?oM)(l + W_,{zq))-\ if 0 < Uf < 8 
-{j3o/j3i){l + Wo{zQ))-^, if9<nf<16 


( 2 ) 


where zq = —/^o and (3i are the first two /3-function coefficients, 
= f^i/Po and A = Aj^ is the conventional MS scheme QCD parameter. Using 
formula the analytical structure of the 2-loop coupling in the complex plane was 
determined imaEi. The motivation for these studies was a need for the development of 
dispersive methods [6-25]. The dispersive approach has been devised to extend properly 
modified perturbation theory calculations towards the low-energy region 1313US]- The 
most simple and elaborated variant of the dispersive approach, the Shirkov-Solovtsov 
Analytic Perturbation Theory (APT), was formulated in refs. [3 and jH] (for a review 
on APT and many original references see refs. PI and [21] )• If should be remarked 
that in the time-like region APT is equivalent to the “contour improved perturbation 
theory” proposed previously in ref. [23 (see also refs. [27-31]). The relation between this 
framework and APT was discussed in refs. II21E2]. More sophisticated nonperturbative 
modifications of the (minimal) analytic QCD model of Shirkov and Solovtsov were also 
presented [18-23]. A generalization of APT to non-integer (fractional) powers of the 
running coupling has also been proposed and successfully used to calculate the three- 
point functions in QCD |23j . 

The 2-loop explicit solution m was soon found to have a more important 
application. In ref. [22] i the running coupling in the k-th order {k > 3) in a MS-like 
renormalization scheme was expanded in powers of the exact two-loop order coupling 
(here and hereof we omit the argument Uf) 


n=l 

On this basis, author of ref. 1331 has proposed a new method for reducing the scheme 
ambiguity for QCD observables. A similar expansion (motivated differently for an 
observable depending on a single scale) was suggested in ref. [21] • Note that the 
analogical expansion but in powers of the approximate (asymptotic) 2-loop coupling 
has previously been introduced in ref. |22j- This expansion was used to construct the 
running coupling with consistent matching conditions at the quark thresholds in the 
3-loop order. However, if the asymptotic 2-loop coupling is used, the coefficients of the 
expansion depend on the scale Q^. The main advantage of dBD is that it allows us to 
write QCD observables, in massless renormalization schemes, as series in powers of the 
renormalization scheme independent quantity, the exact explicit 2-loop coupling ©• 
One could introduce a similar expansion in powers of the 1-loop scheme independent 
coupling as well. However such a series would not be useful, since it could not describe 


t Here = —q^ = —(q°)^ -t ^ and > 0 in the Euclidean domain. 
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the double logarithmic singularities of the higher order coupling. Recently, the series m 
has been used to construct exact explicit expressions for Euclidean and Minkowskian 
observables within APT isainni (see also ref. inn). In practice, the hrst few terms in 
series give the excellent approximations to the coupling even in the infrared region 

[saisE] 

The main purpose of this paper is to present a detailed mathematical investigation 
of the series ©• In Sect. 2 we discuss in some details the singularity structure of the 
MS coupling, in higher orders, in the complex plane. In particular, we determine 
the locations of the Landau singularities of the coupling (at the 3- and 4-loops) as a 
function of rif for the Uf values into the range of validity of the asymptotic freedom of 
QCD. A similar investigation (but for large Uf values when the /3-function has a positive 
hxed point) has previously been undertaken by the authors of refs, [ini Ej using a 
different technique, whose work does not overlap the material in Sect. 2 to a marked 
extent. In Sect. 3 we prove that the series © in the MS-like renormalization schemes 
has a positive radius of convergence to all orders in perturbation theory. In the proof 
we use the methods of the analytical theory of differential equations. In Sect. 4 we 
solve a higher-order RG equation for the running coupling implicitly as a function of 
the 2-loop running coupling. By means of the obtained transcendental equation, we 
determine the analytical structure of the higher order coupling in the complex 2-loop 
coupling plane. As a result, we evaluate analytically the radii of convergence of series © 
at 3- and 4-loops as a function of Uf. In Sect. 5 we determine the convergence region of 
the series solution with respect to the momentum squared variable Q^. For sufficiently 
large Uf values (rif > 14 in the MS scheme), we hnd that the series converges for all 
> 0. For the lower Uf values, we determine the minimal value above which 
the series converges. We compare this scale, at rif = 3, with the infrared boundary 
of perturbative QCD estimated within two different nonperturbative frameworks. In 
Sect. 6 we study properties of the dispersive images of the series solution ©, the non¬ 
power series determined in the sense of the QCD Analytic Perturbation Theory approach 
of Shirkov and Solovtsov, both in the space- and time-like regions. Our conclusions are 
given in Sect. 7. In the Appendix we collect some relevant formulas which we need in 
our computations. 

2. The Analytic Structure of the Coupling to Higher Orders 

In this section, we will determine the location of the Landau singularities of the coupling 
in the complex plane at the 3-loop and 4-loop orders. As we shall see, there is a 
close relation between these singularities and the convergence properties of the series 
©. For large Uf values $, we will reproduce part of the results of jTTn ITTj using another 
technique. Let ns hrst give some familiar aspects of the RG equation in the way we prefer 

f Note that QCD as a realistic theory appears only for 6 flavours and below. However, there are 
different theoretical motivations to consider the multi-flavour theory for 6 < ny < 16 (see, for example, 

[23EH1ISSI) 


A novel series solution 


4 


to formulate it. It is more convenient to introduce the dimensionless variable u = Q^/A'^ 
and a modified running coupling, a{u) = Poas^Q"^), satisfying the RG equation 


u- 


da{u) 

du 


k-l 




-,n+2, 


u 


n=0 


( 4 ) 


where Po) and bn = Pn/Po^^ (for our notations see the Appendix). 

The MS scheme values of the hrst three coefficients &i, 2,3 are listed in Table A.l. 
With this normalization of the coupling relevant perturbative formulas in higher order 
applications of renormalization group become simple 0 cni I2Z|. Suppose that the 
solution to 0 a{u) has a singularity at some hnite point, u = ui^, in the complex u- 
plane, i.e. a{u) —*• cx) as n —>• ml- If follows then from the differential equation Q 
that 


^ /{u — as n —>• 

However, to conhrm the existence of the singularities and to determine their positions, 
a detailed investigation is required in each hnite order of perturbation theory. 

Let us integrate Eq. © for sufficiently large real positive values of u = exp(f) in 
the neighborhood of a real point Uq = exp (to) 

t = Inn = T^^^(a) where T^^^(a) = f {P^^\a')}~^da'+ to, (5) 

J ao 

this can also be written 

t = + 6i ln(a) + T^^^(a) + to, (6) 

where is a regular at zero function 


f^^\a) = / G^'‘\apda' : G^^\a) = 1/P^^\a) + l/a^ - fei/a. 


( 7 ) 


' ao 


here the integration constant has been redehned: to = to — — 6iln(ao). The 

conventional dehnition of the scale A parameter [lU] leads us to the condition to = 
_j’(fc)(0) §. With this choice Eq. (0) reads 


t = a ^ + 6i ln(a) + / G^^\a')da'. 

Jo 

We could write in place of © another but related formula jH 

t = a~^ — bi ln(6i + a“^) + f 

Jo 

where 


( 8 ) 


= (1 +6ia)-‘ j;*. 


the elementary functions. In the 3-loop case, we can write 
T(3)(a) =a-^ + bi In (a) + Ti^\a) - T® (0), 


\a')da' 


( 9 ) 

/k-l 

\ fc-i 



Y, K<r-\ 

( 10 ) 

\n=0 

) n=2 


m T(^)(a) 

can be expressed 

in terms of 


( 11 ) 


§ the term proportional to l/ln^fw) in the asymptotic expansion of a{u) at large u should be suppressed. 
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where 



, . / I"?'!/ 262 — fe? fbi + 2b2a\ 

—O.bbi ln(P® (a)) + arctan if 0 < n/ < 5 

ln(a - ai) ln(a 2 - a) 

7-777-;-7 “ 7-777-;-7 it 6 < Tif < 1 ' 

(oi — a 2 )(l + bitti) (oi — a 2 )(l + ^102) 


( 12 ) 


F(3)(a) = + b^a + 1, A® = 462 - bj, and ai ,2 = (-&i ± V^^)/{2b2). In the MS 

scheme A® > 0 (< 0) if 0 < n/ < 5 (6 < n/ < 16) (see Tables 1 and 2). 

Let us specify the locations of the roots of the algebraic equation 


k-l 

FW(a) = = 0, (13) 

n=0 

the non-trivial zeros of the (3 function, for different values of nf. For 0 < nf < 7, in the 
4-loop order in the MS scheme equation m has one negative real root oi < 0 and a 
pair of complex conjugate roots 02 = 03 (see Table 3). We will assume that 7 ^( 02 ) < 0. 
Then Eq. (jHj) can be written as 

t = a-^ + bi ln(a) + T^^^V) “ ^i^^(O), (14) 

where is a regular at zero function 

= —b^^ {El ln(a — oi) -|- 3 fJ(F 2 ) ln[(a — 02 ) (a — 03 )] 

+29(i!.) arctan } . (15) 

with 


Fj — {flj {cLi 0>j){p>i ®A:)} t ^ 7^ J 7^ ^ 

and {i,j,k) is a cyclic permutation of (1,2,3). Equation (fT^ was derived for real 
positive values of a. It may be continued analytically in the complex a-plane choosing 
the relevant branches for each elementary function on the right hand side. For unphysical 
values 8 < ri/ < 16, in the 4-loop case Eq. has three real roots: oi < 0, 0 < 02 < 03 
(see Table 4). Let a be positive lying in the interval 0 < a < 02 . Formula (jH)) may now 
be rewritten 

3 

t = a~^ + bi ln(a) — Fj ln(a 7 ^(cij — a)). (16) 

i=l 

For complex values of a, the analytical continuation of Eq. (unD can be easily performed 
assuming that each logarithm on the right of have its principal value. 

We have to determine the analytical continuation of the coupling starting from the 
implicit solution ©• This implies a preliminary study of the analytical properties of the 
inverse function t = T^^\a) in the complex coupling plane. The one-valued branch of 
this function may be dehned in the cut complex a-plane choosing the cuts carefully. The 
physical branch may be determined from the requirement that the branch yields a real 
positive t for real positive and sufficiently small values of a. It is convenient to determine 
the analytical continuation of (a) starting directly from the integral representation 










A novel series solution 


6 


Im(a) Im(a) 


Cli 

las 

Re (a) ai 

02 a3Re(a 


y a2 



(a) 0 < Uf <7 (b) 8 < Uf < 16 


Figure 1. The singularity structure of the function t = T^^\a). Two different 
situations are shown. Branch cuts are represented by bold lines and branch points 
by the blobs. 


(jni). The integral there should be regarded as a line integral in the complex a-plane. The 
line must be deformed to avoid singularities of the integrand. It is seen from (jH)) that 
the function has a simple pole as well as a logarithmic branch point at a = 0. In 

addition, there are logarithmic singularities at the roots of the algebraic equation m- 
For 0 < rif < 5, in the 3-loop case (in the MS scheme) Eq. (I13jl has a pair of complex 
conjugate roots, while it has two real (positive and negative) roots for 6 < Uf < 16 
(see Tables 1 and 2). The 4-loop case has already been discussed above. We must 
make a branch cut along the negative a-axis {a : —oo < a < 0 } corresponding to the 
logarithmic branch point at zero. In the cases, where there is a pair of complex conjugate 
branch points (say 02,3 in the 4-loop case) we must choose additional branch cuts. One 
possibility is to choose the cuts parallel to the imaginary axis (see Fig. 1(a)), so that 
these branch points are connected by the cut running through the infinity. In the 4-loop 
case, we choose the cuts {a : —00 < 9(a) < 9 ( 02 ), 9 ( 03 ) < 9(a) < cx),3fJ(a) = 3fJ(a3)} 
(see Fig. 1(a)). The analytical continuation of in the cut complex a-plane will be 

determined uniquely if we require that (a) is real for a real positive and sufficiently 
small values of a. Note that the above considered choice for the cuts is not unique. We 
could, for example, choose the cuts running along straight lines connecting the complex 
conjugate branch points to the origin. Nevertheless, former possibility (which we accept 
in this paper) seems to be preferable: with this choice f as a function of the phase of a 
will be continuous in the neighborhood of a = 0 with the exception of the cut running 
along the negative a-axis. 

Consider now the theoretical cases with only real roots. Let ai be the negative 
root, and 02 be the positive one (in the 4-loop case 02 is the smallest positive root). 
Then the branch cuts may be chosen along the real intervals {a:—cx)<a<0} and 
{a : a 2 < a < cx)} (see Fig. 1(b)). To determine the physical branch, we require that 
T^^\a) is real in the real interval a G (0, a 2 ). We may now analyze the singularity 
structure of the running coupling a = a{t) = a{u) in the complex f-plane, and hence in 
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the complex w-plane too. Evidently, the singular points are determined by the limiting 
values of the function as a tends to inhnity ||. In general, the limiting values may 

depend on the way along which a tends to inhnity. Consider, for example, the 3-loop 
case for 6 < rif < 16. We start from formula (HH). Let a be a point in the complex 
plane a = |a|exp(*(5), where S is the phase of a. The analytical continuation to this 
point gives 


n tr.r'i _ r 

= exp (—Z(5)|a|“^ — 0.56i ln(|P'^^^(a)||a|“^)- 

V —A® 


— ®l) 


—ai(a 2 — a) 

+ i{bi5 - 0.56i(5i + 52) - (0.56? - 62)(5i - 


(17) 


where = 620 ^ + Wa -|- 1 (62 < 0 ) and 6 i {i = 1 , 2 ) denote the increments of the 

arguments of (a — ai) and (02 — a) as a goes from 0 to the point a along a contour T: 
5i = Ar arg (a — Oi) and 62 = Ap arg(a 2 — a). We can now calculate limiting values of 
T(3)(a) as a tends to the sides of the branch cuts along the real axis. We hnd, at the 
sides. 


5 = 0, 5i = 0, 52 = =F7r, if Re(a) > 02 , Im(a) = ±e 

5 = ±7r, 5i = ±7r, 52 = 0, if Re(a) < oi, Im(a) = ±e 

where e —> O’*" is assumed. Using Eq. m with Eq. (HHD, we may easily calculate the 
limits of when a goes to inhnity along the upper or lower side of the right 

(left) cut. One may conhrm that the result will be the same regardless of the branch 

cut chosen. It depends only on the side of the cut: = lima^oo (a ± te) = 

hma^_oo ± te). Thus we hnd, for 6 < uj < 16, 

= -0.56iln|62| + In [ 02 / 01 1 ±i (o.5bi + 

Consider now the 4-Ioop case for 8 < uj < 16. Using arguments similar those employed 
in the 3-loop case, we hnd the singular points 

3 

t^±\n{) = 63^ ^Efcln lafcl ± ^(^2 + Es)^. (20) 

k=l 

Now it is sufficient to show that the obtained limits do not depend on the special choice 
of the directions in the complex a-plane. To see this, let us take the contour integral 
in Eq. m along a closed contour chosen as follows. Let the contour consists of two 
diherent curves with common starting point at zero. Let both curves lie in the upper 
(lower) complex plane and are connected by the arc of a circle with centre at zero and 
radius R. The integrand has no singularity inside the contour, and the value of the 
integral is therefore zero (the Cauchy’s theorem). Consider the limit when the radius 
of the circle tends to inhnity. Then the integral along the arc tends to 0. So that the 
integrals along the two different curves tend to the same limit. Thus the result stated 
follows. 

It is important to determine whether or not the singular points t^±\n{) are located 
inside the strip —tt < 7s(t) < tt. The strip is image of the hrst (physical) sheet of the 


II we assume that the coupling does not have singular points where it takes finite values. 
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complex plane under t = ln(Q^/A^). Depending on the value of rif there are two cases 
to consider. The hrst case is that the points lie inside the strip so that the unphysical 
Landau singularities appear in the hrst sheet. Then the running coupling is not causal, 
and thus perturbation theory is incomplete: the non-perturbative contributions are 
required to remove the unphysical singularities p El ca EE]. This case corresponds 
to real-world QCD, where rif < 6. In the second case, the singular points may arise 
beyond the strip. So that there are not real or complex singularities on the hrst sheet 
of the momentum squared variable, and thus perturbation theory is consistent with 
causality. The singularities still may present only on the time-like axis < 0. The 
second possibility may be realized for sufficiently large Uf values. The value of Uf above 
which the causal analytical structure of the coupling is restored can be found from the 
equation 

= ±7r. (21) 

With the 3-loop formula (ITT?|l . we hnd the solution to (pT|l 8.460, and with the 

4-loop formula (HB ^ 8.455. Thus, the 3- and 4-loop MS scheme results almost 
coincide. We remark that the 3-loop estimation was obtained previously in EH]- The 
2-loop condition for causality of the coupling can be found in p. In our notation it 
reads &i(nf) < —1, this gives for the lower boundary of the causal region the value 
^ 9.68. 

Note that for Uf > the /d-function has a positive infrared stable hxed point (see 
Tables 2 and 4). So that the running coupling is trapped in the range between 0 (the 
ultraviolet hxed point) and the infrared hxed point at all energies. The fact that QCD 
in perturbation theory for sufficiently large Uf values has an infrared hxed point has 
long been discussed EZ|. Of particular interest is the case when the hxed value of the 
coupling is sufficiently small. Then, presumably, the theory may be reliable described 
within perturbation theory for all energies including infrared region. The corresponding 
interval of Uf values is called as a conformal window P EEl EH EH!- It is believed 
that there is a phase transition in QCD with respect to Uf inside the range of validity 
of asymptotic freedom 0 < Uf < 16. For small values of Uf below the critical point 
(uf < < 16) the theory is dehned via the conhning phase. Above this point, there 

is a conformal window < Uf < 16, where the theory is dehned via the non-Abelian 
Coulomb phase with no color conhnement and dynamical chiral symmetry breaking. 
One possible way to determine the critical point is to use the Oehme-Zimmermann 
criterion for the gluon conhnement, the superconvergence rule for the transverse gluon 
propagator [^. This gives the value A"" = 13Ac/4 (=9.75 for W = 3 colours). Other 
possibility is to apply arguments of dynamical chiral symmetry breaking P EH 1121 IHj ■ 
This gives slightly large value A^" ~ 4A),. Assuming the value for as predicted 
by the superconvergence rule, the authors of EHl EZI have given arguments that the 
perturbative running coupling inside the conformal window (and beyond the 1-loop 
approximation) is always causal, i.e. < A"". Note that the infrared hxed value of 
the coupling inside the window is not large: it coincides with the root 02 (in the 3-loop 
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case, for Uf = 10, 02 ~ 0.26). 

Consider now the cases where Eq. (USD has complex roots. In the MS scheme this 
takes place at 3-loops if 0 < rif < 5 and at 4-loops if 0 < rif < 7. The corresponding cuts 
in the a-plane are chosen as shown in Fig. 1(a). We hrst calculate the limit of T^^\a) 
as a tends to inhnity along a line going to inhnity in the right half-plane 3fJ(a) > 3fJ(a2). 
Evidently, the result will not depend on the particular choice of the direction as far as 
the line belongs to the right-half plane 3fJ(a) > 3fJ(a2). Choosing the path along the 
positive semi-axis and using the 3- and 4-loop formulas (USD and (USD, we calculate the 
limits, = lim^^oo(a). 


^!hpK) = -0.56i ln( 62 ) + 
= ^3 ' (- 255 (^ 2 ) 


262 - bj 
Vaw 


/ TT 

- arctan 

V2 


^O.Stt -|- arctan 


^ia2) 

|G(a2) 


+Ei In |ai| -h 23?(E2) In |a 2 |), 



( 22 ) 

(23) 


here the subscript “rhp” shows that the limits are calculated along the way going to 
inhnity through the right half plane 3?(a) > 3 ^( 02 ). Let us now calculate the limits 
of when a tends to inhnity through the left half plane 3fJ(a) < 3 ^( 02 ). We 

may take without loss of generality the ways along the sides of the cut running on 
the negative semi-axis. The limiting values of from above and below the cut, 

T^\a) = lim£^o+ ±?e), may be determined by the analytical continuation of the 

right hand sides of Eqs. (ED and ED- The singularities are then determined by the 
limiting values = lima^_oo (a). This gives 

t[JJ,±(nf) = - ^(2^2 - fe?)/VA® ± *7r6i, (24) 

= ^Ap(%) + 27r3(E2)/63 ± - Ei/bs). (25) 

Equivalently, by the analytical continuation of m for negative values of a we obtain the 
useful formula 

/*—oo 

=-bi\nbi+p.v. / g''’‘\a')da'± iTr{bi + Tes[g^^\a), (26) 

Jo 

where ^ denote the real negative root of (fT!?|l (which presents only in the 4-th order 
case) and the integral here is considered in the “principal value sense” (p.v.). It is seen 
from the Tables 1 and 3 that in these cases the Landau singularities present in the hrst 
sheet of the Q^-plane. 


3. The Proof of the Convergence of the Series 

Inserting series dSD into the RG equation dUD, we recursively determine the coefficients 
in terms of C 2 (ci = 1) and the /3-function coefficients. However, the coefficient 

(k) 

C 2 ^ still remains undetermined. This rehects the arbitrariness in the dehnition of the 

(k) 

A-parameter. With the conventional dehnition of the parameter, we hnd that C 2 = 0. 
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Table 1. The 3-loop MS quantities: the complex zeros of the /3-function 01^27 and the 
(3) (3) 

singularity locations and in the t-plane as a function of nf for 0 < nf < 5. 


nf 

01.2 

‘'rhp 

hhv± 

0 

-0.393 ± 0.8821 

0.844 

-1.539±2.648z 

1 

-0.399 ± 0.8921 

0.839 

-1.506± 2.628z 

2 

-0.415 ±0.916z 

0.830 

-1.433 ±2.578z 

3 

-0.447 ±0.965z 

0.810 

-1.293 ±2.482z 

4 

-0.526 ± 1.071Z 

0.766 

-1.026 ±2.322z 

5 

-0.819 ± 1.349Z 

0.650 

-0.423±2.067z 


Table 2. Same as in Table 2, but for 6 < nf < 16. 


nf 

Oi 

02 


6 

-1.49 

7.09 

0.173 t0.077z 

7 

-0.877 

1.24 

-0.159 T l-05z 

8 

-0.651 

0.660 

-0.019 t2.36z 

9 

-0.509 

0.409 

0.620 t4.26z 

10 

-0.403 

0.264 

2.17t 7.21z 

12 

-0.245 

0.104 

13.6t21.3z 

14 

-0.125 

0.028 

98.7 =F 89.9z 

16 

-0.023 

0.001 

6216t2852z 


Table 3. The 4-loop MS quantities for 0 < nf < 7. 


nf 

Ol 

02.3 

‘'rhp 

''lhp± 

0 

-0.797 

0.130 TO.782z 

1.164 

-1.294 ±0.961z 

1 

-0.794 

0.134 TO.784z 

1.163 

-1.284 ±0.930z 

2 

-0.793 

0.142 TO.792z 

1.158 

-1.256 ±0.871z 

3 

-0.796 

0.159 t0.810z 

1.145 

-1.199 ±0.768z 

4 

-0.802 

0.190 t0.844z 

1.114 

-1.099 ±0.593z 

5 

-0.806 

0.259 t0.906z 

1.035 

-0.940 ±0.295z 

6 

-0.795 

0.444 T 1.012Z 

0.821 

-0.707 t0.217z 

7 

-0.741 

1.124 t0.973z 

-0.047 

-0.423 T 1-074Z 


Table 4. Same as in Table 4, but for 8 < nf < 16. 


nf 

Ol 

02 

03 


8 

-0.623 

0.699 

6.476 

-0.093 t2.37z 

9 

-0.482 

0.426 

4.810 

0.522 t4.28z 

10 

-0.362 

0.281 

1.937 

1.95t7.32z 

12 

-0.193 

0.112 

0.548 

13.1t21.9z 

14 

-0.085 

0.030 

0.196 

98.1t91.8z 

16 

-0.014 

0.001 

0.031 

6219t 2858z 
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This follows from Eqs. m and 0 if we use the asymptotic expansion for the Lambert-hh 
function (see pp. 22-23 in |3]) 

Wk (z) = Li-L 2 + L 2 /L 1 + L 2(-2 + L 2 )/ 2 LI + L 2(6 - 9 L 2 + 2LI)/QL\ + 0 {{L 2 /Uf) (27) 

where for the branch W-i{z) for real negative z {z ^ 0~ as —>■ cx)) we must put 

Li = ln(— 2 ;) and L 2 = ln(—ln(—z)). Several hrst coefficients calculated in the 4-loop 
case are: = 1, = 0, = /d2//9o, = l3‘i/2j3Q, 

(4) _ 5^ ^ <dl/?3 (4) __1 1 , ^ /^2<d3 

~3^l “ 12 PI "^12 PI ^ PI ’•••• 

Inserting these values of the coefficients into series ® and using (EZD, one may readily 
reproduce the conventional asymptotic representation for the coupling (see, for example, 


a 


asy 


(Q^) = 


Pi In L 


+ 


AL HI n 




InL 


■’*1 


+ 0 


In^L 

L4 


(28) 


where L = \n{Q '^/> 1. 

Let us change the variable according to ^ 6 = PQaf\Q‘^) and introduce the 
new function w{6) = w^^\9) = {Q"^)/af\Q‘^) — 1. The RG equation (jA.lfl may be 
rewritten as 




(29) 


where 


/“’(». •»)= i! '>"""(1+•»)"-(•»+1). (30) 

with bn = Pn/Po~^^- The function of two variables f^^\w,9) has the Taylor expansion 

00 

/<‘)(«,, 9 )= ^ ( 31 ) 

m,n=0 

with ? 7 q^q = 0, 77(^0 = 1 Vo^i = 0- th® 4-loop case, the expansion is 

9) = w + bi9w + 626 *^ + w‘^ + (63 - bib2)9^ + - (32) 


We may now use the analytical theory of differential equations jini EE) to investigate 
Eq. (pHl . With the initial condition ta(0) = 0, this equation has a singularity: for 
0 = 0 and w = 0 the ratio f^^\9,w)/9 is undehned. Nevertheless, in the special case 
where = 0, = 1 and = 0, Eq. (1^ may still have an analytic solution 

satisfying the initial condition t(;( 0 ) = 0 (see e.g. and IIEI)- For the sake of clarity, 
the following discussion is quite detailed. The expansion (ED) converges in the domain 
D = {0 < |ta| < ri, 0 < | 0 | < r 2 }, where ri and r 2 are some positive numbers chosen in 
the range {ri,r 2 : ri < 00 , r 2 < l/| 6 i|}. It follows then from the classical theory that 
there exists a positive number such that \ f^^\w,9)\ < for {w,9) C D, and 
the coefficients rjm^n satisfy the inequalities jUl EOl 

I’ffnI < 4t’„. Where 


(33) 
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Under these conditions, we will show that there exists a regular solution to Eq. (1^ 


u>(«) = «>'*>(») = 


(34) 


n=2 


where Cn^ = with Cn'^ being the coefficients in the original series (jSI). We 

recall that according to our choice = 0. Inserting expansions m and 

dSl into Eq. (i2ni) we recursively determine the coefficients 


•-2 — 'lo,2y ^'-3 — 'll,1^2 “r '/o,3 ) ' ' ' 


(35) 


(36) 


(37) 


Consider now the auxiliary function w = w{9) satisfying the equation 

w = 9) = - - ^ ^ U + — + — 

[1 - w/ri)[l - 9/r2) V ^2 

it has the Taylor expansion 

m=0,n=0 

with the coefficients dehned in Equation (ESI) has a series solution 

OO 

w{9) =w^^\9) = 

n=2 

Inserting the expansions (EH) and (EHI) into Eq. (ESI) we hnd the recurrence formulas 

7f=dU?’(39) 


(38) 


Let us compare Eqs. (ESI) with Eqs. ESD- Making use of Eq. ESD, we obtain the 
inequalities \cn\ < for n = 2,3.... It follows then from the comparison test 

that the series El is absolutely convergent in the disk of convergence of the series (IHHjl . 
Evidently, the series ESD has a positive radius of convergence. The radius is equal to 
the modulus of the nearest to the origin singularity 9^r of the function w = w{9). The 
value w can be solved explicitly from the quadratic equation ESD- The singularities of 
the majorant function w{9) are therefore located at the zeros of the discriminant of the 
quadratic equation. Hence we hnd the singular points 

0f)=r2, 9f^ = -rir2/M^^\ 9^} = /2 - 1 ± \/m^ + 8/), (40) 

where = 1 + To obtain the best possible estimation, we have to maximize 

|9S’| with respect to ri and r 2 . The quantity = maxu, ^ |/(^)(tc, 6*)| may be 

determined according to the maximum modulus principle. The modulus \f^^\w,9)\ 
takes its maximum on the circles |tc| = ri and l^l = r 2 . We hnd that the maximum 
of |/^*^)(w,0)| is attained for real positive values of w and 9. Furthermore, 9^^ = 9^^\ 
We choose the values ri = 0.25 and r 2 = 0.42 in the 3-loop case, while ri = 0.21 
and r 2 = 0.32 in the 4-loop case. Using “Maple 7” E3, determine numerically 
the maximal values of \f^^\w,9)\ on these circles. For rij = 3, we have found that 
0.695 and Ri 0.596. Computing the modulus of the numbers (I40j) and 
comparing them, we determine the radii of convergence of the majorant series (EHD- 
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Pmaj = ~ 0.045 and = 16 * 3 “^^| 0.033. Thus we have found that the radii of 

convergence of the original series (0 in the MS scheme in the 3- and 4-loop orders, at 
Hf = 3, are bounded below as > p^aj/l^o ~ 0.06 and > p^^j/Po ~ 0.05. As 
we shall see, the actual values of p^^l are significantly larger than the obtained lower 
bounds. 

We remark that the above proof of the convergence of the series holds in all MS-like 
(massless) renormalization schemes, since in the proof we have not used specihc values 
of the /^-function coefficients and the condition C 2 = 0 is common for all these schemes. 


4. Determination of the Radius of Convergence of the Series 


By a change of variable ^ 9 = a^'^\u) {u = Q^/A^) Eq. (j3|) can be rewritten 

da _ 


(41) 


de e^ + hiO^ ’ 

where a = F^^\9) = {u) = f3QQ^\Q‘^). In the following we will sometimes, but not 

always, omit the superscript “(/c)” referring to the order of perturbation theory. In the 
preceding section, we have shown that the series (ED or equivalently the series 


a = F{e) = Y,cne^^ (Cn = A 


—n+1 ^ 


(42) 


n=l 


has a positive convergence radius. It is possible then to dehne the inverse function 
9 = F~^{a), which can be expanded in powers of a 


9 = F 


-1 


00 


(fl) ^ ^ dfiCL 


(43) 


n=l 


By using arguments similar to those employed in Sect. 3, one can verify that the series 
(gSD also has a hnite radius of convergence. Under this condition, we may apply the 
classical method for estimating the convergence radius of series (see the book E] PP- 
146-148). The main argument is that the function a = F{9) must have at least one 
singular point on the circle of convergence of the series (14211 . There are two possible 
cases. First, suppose that 9o be a hnite singularity of F{9), where the function takes a 
hnite value, oq = F{9o) < cx), while its derivative does not exist. In terms of the inverse 
function 9 = F~^{a) these conditions read 

dF-Ua 


da 


= 0 . 


a=aQ 


Using Eq. dUD at 9 = 9o, we may rewrite Eq. in the form 


d9 

da 


91(1+ 

n+2 


n 


= 0 , 


(44) 


(45) 


a=ao 

for a hnite Oq (which is not a root of = 0 ) has only two solutions: 

6^0 = 0 and 


9q = —1/bi (= —81/64 for Uf = 3). 


(46) 
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Figure 2. The two real branches of W(z). The solid line corresponds to Wo(z) and 
the dotted line in the range W < —1 refers to W-i{z). 


The solution = 0 must be rejected, since at 9 = 0 Eq. iiD does not hold because 
of the initial condition F{6)/6 —> 1 as 0 — 0 (see (jl2))). Secondly, suppose that there 
exists a curve C going to inhnity in the domain of analyticity oi 9 = F~^{a) such that 
F~^{a) 9s < oo as a tends to inhnity along this curve. Then, F{9) has a singularity 
at 6^ = 9s. 

First, we shall consider Eq. m in the 3- and 4-loop orders for 0 < nj < 5 and 
0 < n,/ < 7 respectively. Let us integrate Eq. PT|1 in the real range {9,a : 9 > 0,a > 0}. 
We write the result in the symmetrical form 

1/0 — 6i ln(6i + 1/6*) = 1/a — 6i ln(6i + 1/a) + f g^^\a)da', (47) 

Jo 

where the function g^^\a) is dehned by m, and we have determined the integration 
constant according to the previous choice (see Sect. 2). Equation (ITfl) may be continued 
for complex values of a and 9. Then the integral with respect to a should be regarded 
as a line integral in the complex a-plane. The contour connecting the origin to a must 
avoid singular points of the integrand. 

Note that all the coefficients in series (m are real, so are the coefficients in the 
inverse series fj43j) . Therefore, 0 as a function of a must be real for real values of a 
provided that a is any point inside the circle of convergence of the inverse series. It is 
seen from Eq. iH) that there exists a real neighborhood of the origin 0 = 0, where 
the function F{9) is real and strictly increasing. Thus, in the 3-loop case, the derivative 
F\9) > 0 if 0 > —1/01, provided that 0 < ri/ < 5 (6i > 0). In the 4-loop case, the same 
is true if 0 > max(—l/6i, a|^^^), where is the real negative root of the equation (IT!?|l . 

= 3) ~ —0.796). From this with the initial condition F“^(0) = 0, it follows 
that there exists a sufficiently small real interval including the origin, where the function 
0 = F~^{a) is real and strictly increasing. Inside this interval F ^(a) > 0 (F ^(a) < 0) 
if a > 0 (a < 0). Fortunately, we may solve the transcendental equation (jTfjl for 0 
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explicitly as a function of a in terms of the Lambert-W function 

e = F-\a) = - 6 r'(l + Wn(z))-\ (48) 

where z = ({a) = —(e 6 i)“^ exp(—T(a)/ 6 i), and 

T(a) = — 6 i ln( 6 i + a“^) + f g{a')da'. (49) 

Jo 

It follows from the above discussion that in the region a > 0 inside the convergence disc 
of the series (where 0 > 0) the required branch in Eq. is so that 

e = F-i(^)(a) = -b^\l + (50) 

Formula dSOD may be continued beyond the convergence circle on the positive a-semi- 
axis. It follows from (im that the function z = ({a) is negative and monotonically 
decreasing in the inhnite interval a G (0, cxo). For the considered values of n/ (with 
the exception of Uf = 7 case at 4-loops) ({a) takes values in the range (—e“^,0“) for 
a G (0, oo). The branch W_i{z) is real and negative with W_i{z) G (—cxd, — 1) for 
z G (— e“^,0“) (see Fig. 2). Therefore, the function 6 = F~^{a) determined by (jKn|l 
is real and positive in the entire positive a-semi-axis. However, in the 4-loop case at 
Uf = 7, we hnd that 2 ; G (— e“^, 0“) only inside the interval a G (0, Ob) where Ob ~ 1.003, 
and 2 ; < —e~^ if a > a^. Since W_i{z) has a branch point at 2 ; = —e“^, the function 
F“^*^^^(a) at rif = 7 will have a branch point at a = Ob. The corresponding branch cut 
may be chosen along the positive interval {a : Ob < a < cx)}. Formula PI|) still holds 
on the upper side of the cut: on the sides of the cut the function at Uf = 7 

takes complex conjugate values. Using dSDD, we calculate the limit of F ^(a) as a —> cxd 
along the positive a-semi-axis. So we hnd a singularity of the function a = F{6) 

eS = + fr_i(ci‘>(oo))-‘ = 0^1(4;,), (51) 

here we have used the formula 

C(")(oo)= lim c(U(a) = -(5,e)-'exp(-T«(cx))/6i) = -( 6 ie)-i(n;jj,)-i/'b (52) 

a^+oo ^ 

where = exp (t^hp) being the Landau singularity located on the positive a-semi-axis 
(see the 3-and 4-loop formulas and (I2S1))- The last equality in (E2D follows from 
Eq. (jni). In the 4-loop order, at n/ = 7, we hnd a pair of complex conjugate singular 
points, = lima ^+00 0 (a ± ? 0 ), where 9s.i+ is determined by (|KT|l . 

Formally we may continue Eq. (jUD for negative real values of the variables in 
the region {a, 9 : —l/5i < a < 0,—1/bi < 6 ^ < 0}. Assuming that each logarithm in 
Eq. (I47|l have its principal value, we obtain in this region the equation 

1/0 — 61 ln(l/| 6 '| — 61 ) = 1/a — ln(l/|a| — 61 )-|- f g{a')da'. (53) 

Jo 

Note that the right-hand side of Eq. dSSl) is in fact regular at a = — 1 / 61 : the logarithmic 
singularities of the last two terms are cancelled in the sum. For negative values of a, the 
path of integration of the integral on the right of (jiSl) goes along negative a-axis, but 
avoids the poles of g{a) by small semi-circles above or below the axis. Equation (j5^ 
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Table 5. Locations of the extra singularities in the a-plane at the 3-and 4-loop orders. 


nf 


l4±l 

W 

l«b±l 

0 

0.038 ±0.732* 

0.733 

0.156 ±0.600* 

0.616 

1 

0.038 ±0.741* 

0.742 

0.159 ±0.600* 

0.620 

2 

0.039 ±0.762* 

0.763 

0.166 ±0.609* 

0.632 

3 

0.042 ±0.806* 

0.807 

0.179 ±0.629* 

0.654 

4 

0.046 ±0.902* 

0.904 

0.207 ±0.668* 

0.700 

5 

0.060 ± 1.186* 

1.187 

0.275 ±0.744* 

0.793 

6 



0.498 ±0.888* 

1.018 

7 



1.003 

1.003 


has exactly one real negative solution for 6 inside the interval a E (a, 0), where a = —oo 
in the 3-loop order and it is the hnite negative root of m in the 4-loop order. This 
solution is determined in terms of the branch Wq{z) 

9 = = -b^\l + Wo{C^^\a)))-\ (54) 

where = (e6i)“^ exp(—T(^)(a)/6i) and 

= a~^ — bi\n{—a~^ — bi) + ( g^^\a')da'. (55) 

Jo 

It is instructive to check that our choice for the branches on the real a-axis in fact 
follows from the analytical continuation. To see this, let us expand expressions dsni) and 
dsi as a —>• 0+ and a —>• 0 respectively. We must use expansion (jTfjl for W-i{z) as 
2 ; —> 0“. The same formula, but with Li = In 5 and L 2 = In In 5 holds for Wo{z) 
as 5 —> 00 (see formula (4.19) in ref. |^). One may verify that both expansions 
reproduce the same convergent power series (03). So that Eqs. (EOD and (P|) represent 
the same analytical function in two different regions. Let us now discuss the analytical 
structure of the function 6 = F~^{a) starting from Eq. 0H|) . In general, F ^(a) may 
have singularities at the same points where T{a) is singular. Nevertheless, as we have 
shown, F~^{a) is regular at a = 0, while T{a) is singular there. Furthermore, F~^{a) 
may have additional singularities ai,± arising due to the common branch point of ITo(^) 
and W±i{z) at 2 ; = —e~^. To determine locations of these singularities we numerically 
solve the equation 

^ = (W(a) = -e-^ (56) 

at the 3- and 4 loop orders. The approximate locations of these singularities for different 
rif values are given in Table 5. Note that not all roots of Eq. dsni) give rise to the 
singularities of F~^{a) on the hrst sheet of the Riemannian surface. Thus in the 4-loop 
case, for Uf = 7, Eq. (jsni) has the roots ab± ~ 0.16 ± 0.40*. But one may verify that 
these points are not singular on the hrst sheet. To make the function 6 = F~^{a) 
single valued, we must draw cuts in the complex a-plane taking into account the branch 
points of T{a) and those of W function. We draw cuts in the complex a-plane attached 
to the complex conjugate branch points (say 02,3 in the 4-loop case, the roots of 
and running to inhnity parallel to the imaginary a axis: {a : 3fJ(a) = 3fJ(a2,3), —00 < 
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Q{a) < '^{ 02 ), < Q'(a) < cx)}. In the 4-loop case, we must draw an extra cut 

attached to the real branch point at Oi < 0 and running along the negative semi-axis 
{a : —00 < a < ai}. We also make the branch cuts running parallel to the imaginary 
axis {a : 3?(a) = 3fJ(ab ±),—00 < 9'(a) < Q'(ab-), 9(ab+) < Q'(a) < cx)} attached to the 
branch points at ab±, the roots of Eq. dsni)- However, in the 4-loop case, for rif = 7, 
the branch cut must be chosen along the positive real axis {a : 1.003 < a < cxo}. With 
this choice of the cuts, the function 6 = F~^{a) will be analytic in the disc with centre 
the origin and radius Vc = min{|ab±|, |aj|} (a*, i = 1, 2 ... denote the roots of Eq. (P7?|l h 
For 0 < rif < 5, the points ab± are closest to the origin singularities both in the 3- and 
4-loop orders. Let us dehne a = rexp (zd) and 2 ; = ((a) = \z\ exp (^<h), where 

1^1 = (e 6 i)“^ exp(—3fJT(a)/6i) and ^ = tt — 7sT{a)/bi. (57) 

Let a describes the semi-circle of radius r < lying in the upper half plane (0 < 5 < tt) 
starting from the positive semi-axis. Then the image under z = ({a) describes a curve in 
the 2 ;-plane. This curve intersects the real negative z-semi-axis for two or more times at 
different points. The number of the intersections depends on the value of r: it increases 
when r decreases. At the intersections the boundary of the branch of W is reached, 

so that the branch of W-function must be changed when 2 ; passes these points. To 

dehne the analytical continuation along the semi-circle, we demand that the function 
9 = F~^{a) = F~^{({a)) will be continuous as a function of the phase of a. This will be 
achieved, if we use the rules of counter-clockwise continuity | 3 ] for selecting the branches 
of W when the curve crosses the branch cut. These rules are given by 

kE_i(x -I- ?0) = Wi{x — -?0) if — e~^ < x < 0 

kLi(x — ?0) = Ho(x-t-*0) if — 00 < x <—e~^ (58) 

Wn{x + ? 0 ) = HA+i(x — ? 0 ) if — cxD < X < 0 and n > 1 . 

We start at a = r < Tc on the positive semi-axis with the branch W-i{z) and pass the 
semi-circle {5 : 0 < 5 < tt} selecting the relevant branches according to the rules (IKH|l . 
Using “Maple 7” Ba, we have plotted the function z = ((a) along the semi-circles for 
various values of r in the interval 0 < r < Tc. In this way, we have determined the 
variations of the phase <I> = arg( 2 ;) along the semi-circles. Then we have conhrmed that 
the analytical continuation on the negative interval —< a < 0 , with the rules dSHl), 
really leads to the branch Wo( 2 ;). 

Having the analytical structure of F~^(a) established, we can construct explicit 
expressions for F~^(a) in the entire cut complex a-plane. This enables us to calculate 
the limits of F~^(a) as a tends to inhnity along different directions in the complex plane 
and determine thereby the singularities of the function a = F{6) in the complex 6 '-plane. 
Using the arguments based on Cauchy’s theorem (see Sect. 2), one may justify that it 
is sufficient to calculate the limits choosing the directions only along the real a-axis. 

Let us now dehne the analytical continuation along entire negative a-semi-axis. In 
the 3-Ioop case, we may represent dSSD in the form 

T^^^(a) = — 61 ln(||a|“^ — 5i|)-|-p.v. f g^^\a') da (59) 

Jo 
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Table 6. Positions of the singularities in the 6>-plane at the 3-loop order. 


nf 

0 

1 

2 

3 

4 

5 

Oo 

-1.186 

-1.195 

-1.219 

-1.266 

-1.353 

-1.520 

/i(3) 

^s.l 

0.627 

0.635 

0.653 

0.691 

0.776 

1.029 

^(3) 

^s.2 

-0.594 

-0.601 

-0.618 

-0.653 

-0.731 

-0.956 

p 

0.594 

0.601 

0.618 

0.653 

0.731 

0.956 


for all a < 0. Thus is real, and therefore > 0 for all a < 0. Hence 

fTo(C®(a)) > 0 (see Fig. 2), so that F < 0 for all a < 0. It is evident that the 

required branch, in this case, will be IFo(^) on the entire negative a-semi-axis. Making 
a —>■ —oo in Eq. (jSl and using Eq. dsni), we hnd the real singular point 

hm F-h3)(a) = -(6i(l + lEo(C^'n-oo)))-' (60) 

a—OO 

where oo) = {ebi) ^exp{—T^^\ —oo)/bi), and 

fO 

r®(— cx)) = lim T^^^(a) = —6i ln6i — p.v. / g^^\a) da = ,), (61) 

a—»■—OO / 

J —oo 

the last equality in Eq. dnn) follows from Eq. 

In the 4-loop order, the function is real only inside the finite interval (oi, 0) 

of the negative semi-axis, where ai is the root of (|T3|l (oi ~ —0.796 for Uf = 3). The 
function has a logarithmic branch point at a = Oi. The corresponding branch cut may 
be chosen along the inhnite interval (—cx),ai). To continue T^^\a) in the complex a 
plane, we use formula dSSl). The limiting values of this analytic function from above and 
below the left-hand cut are given by 

T].‘^^(a) = 1/a — 6i ln(|6i — l/|a| I)-|-p.v. f g^*\s) ds ± zktt, (62) 

Jo 

where n stands for the residue 


~ (^2 + & 3 ®i){^ 3 (l + biai){ai — a2)(ai — 03)} 


-1 


K = hm (a 

a—>ai 

and Oj, i=1..3, denotes the roots of (HD- Formula (ED enables us to dehne the analytical 
continuation of the function 9 = in the cut complex a plane. In particular, we 

need to calculate the boundary values of F~^^‘^l{a) as Im(a) —0^ along the left-hand 
cut {a : — cx) < a < ai}. It is easy to convince that the required branch of W along the 
sides of the left-hand cut will be Wo{z), provided that \n/bi\ < 1. This condition holds 
only for 0 < rij- < 5 (for example, n/bi 0.690 at nf = 3). Therefore, for 0 < uj < 5, 
formula (ED is valid also along the sides of the cut. But, for 6 < n/ < 7, we have 
\K/bi\ > 1. Then one may check that the relevant branches on the opposite sides of 
the cut should be IT±i. We may now calculate the limit of F~^^^\a) as a approaches 
inhnity going along the upper or lower side of the cut. So we determine the singularities 
of F(^)(6'). We use Eq. (IKip with Eq. if 0 < Hf < 5. However, for Q < Uf < 7, the 
branch Wq in Eq. (ED should be replaced by W±i. Then we find the singular points 

lim F-«‘)(a±,0) = | + if0<n,<5 

“ 1 —fe, ^(1 + (— 00 ))) ^ if 6 < nf < 7 


( 63 ) 
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Table 7. Positions of the singularities in the 0-plane at 4-loops. ni and n 2 denote the 
labels of the branches of the W function used to calculate the singularities 0s.i and 
^s. 2 ± respectively, p is the radius of convergence of the series 


nf 

00 

^s.l 


ni 

0s.2± 

^s.2± 

n2 

p 

0 

-1.186 

0.485 

- 

-1 

-0.545 T0.334^ 

0.639 

0 

0.485 

1 

-1.195 

0.488 

- 

-1 

-0.544 TO.341z 

0.642 

0 

0.488 

2 

-1.219 

0.497 

- 

-1 

-0.546 t 0.354z 

0.650 

0 

0.497 

3 

-1.266 

0.516 

- 

-1 

-0.550 t0.380z 

0.668 

0 

0.516 

4 

-1.353 

0.554 

- 

-1 

-0.552 T0.429^ 

0.699 

0 

0.554 

5 

-1.520 

0.641 

- 

-1 

-0.533 TO.526z 

0.748 

0 

0.641 

6 

-1.885 

0.934 

- 

-1 

-0.394 ± 0.672^ 

0.779 

±1 

0.779 

7 

-3.008 

-0.887 T 1-53L 

1.769 

Tl 

-0.105 ±0.614^ 

0.623 

±1 

0.623 

Table 8. Numerical examination of the series 14211 to 3-th and 4-th orders at nf = 3. 


n 

Pn 

Pn 

n 

~(h 

Pn 

Pn 

10 

1.04 

0.737 

80 

0.692 

0.550 

20 

0.807 

0.632 

90 

0.688 

0.547 

30 

0.752 

0.597 

100 

0.685 

0.544 

40 

0.727 

0.579 

no 

0.682 

0.542 

50 

0.713 

0.568 

120 

0.680 

0.540 

60 

0.703 

0.560 

130 

0.678 

0.538 

70 

0.697 

0.554 

140 

0.677 

0.537 


where C±^(—cx)) = (e6i) ^exp(—cx))/6i), and by Eqs. and (OF)|l 

T^\—oc)= lim T^‘^\a ± tO) = ^{t\hp±) ± iTTK, (64) 

a^—oc 


here the subscript “±” shows that the limits were evaluated keeping the upper (lower) 
side of the cut. Evidently, = o), 2 +. 

In Tables 6 and 7, we tabulate the singularities 9o{ni) = 9^ll{ni) and 9^^}{ni) 

at the 3- and 4-loop orders respectively. In the 3-th order, we see that the singular points 
(S') (s') 

2 are closer to the origin than the points 9o or 9^ {, so that the radius of convergence 

/ o\ 

of the series (m is equal to | 6 'g 2 d. On the contrary, in the 4th order, we hnd that 
= 1 6 ^ 3 j I if 0 < rif < 5, and = 103^2 1 if 6 < rif < 7. We recall that the radius of 
convergence of the original series (jSl) is determined by p^^^ = p^^^Po (p^^^ = 0.965 and 
= 0.720 for nf = 3). 

Looking at the numbers in Tables 6 and 7 one sees that the convergence radii in the 
range 0 < nf < 6 increases as Uf increases. But, for hxed values of n/, they decrease as 
the order of perturbation theory increases. In order to examine the obtained formulas, 
we calculate numerically the coefficients of the series (US) Cn for large values of n. In 
Table 8 , we study the behaviour of the quantity pn^ = 4-loop 

orders at Uf = 3. It is seen from the table that our predictions are in good agreement 
with the expected limiting relation pn'’ —>■ as n —>■ cx). The theoretical predictions 

are p^^^ = 0.653 and p^^'> = 0.516 for Uf = 3 (see Tables 6 and 7). 
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Table 9. The location of the branch points ai_ 2 , Ob and ab± of 6* = F for 

rif = {6 — 16}. 


nf 

ay 

02 

ay, 

ah± 

6 

-1.489 

7.089 

2.418 

-2.258 ± 0.296^ 

7 

-0.877 

1.238 

0.813 

-0.846 iO.lOF 

8 

-0.651 

0.660 

0.541 

-0.552 ±0.006^ 

9 

-0.509 

0.409 

-0.360 

0.439 ± 0.036z 

10 

-0.403 

0.264 

-0.253 

0.266 ±0.034z 

11 

-0.318 

0.169 

-0.177 

0.159 ±0.027z 

12 

-0.245 

0.104 

-0.121 

0.089 ±0.018z 

13 

-0.182 

0.059 

-0.078 

0.043 ± 0.01b 

14 

-0.125 

0.028 

-0.044 

0.015 ±0.006z 

15 

-0.072 

0.010 

-0.020 

0.002 ±0.002z 

16 

-0.023 

0.001 

-0.004 

-0.0013 ±0.0002z 


Next consider the theoretical cases with large rif values where the /3-function has 
non-trivial real zeros. This takes place in the 3-loop case for rif = {6 — 16}. From now 
on we shall conhne ourselves to the 3-loop case. There are now two different cases which 
should be considered separately. For rif = {6 — 8} ( 6 i > 0 and 62 < 0) the equation 
= —e~^ has one real positive root at, (0 < a?, < 02 ) and a pair of complex 
conjugate roots ab± with 3?(ab±) < 0 (see Table 9 ). On the real interval ai < a < ay,, 
the real analytic solution to iZI) is 

j + if 0 < a < Ob . . 

+ Wq{z))~^ if Oi < a < 0 , 

where z = C^^\a) (see Eq. ii). To determine uniquely the analytical continuation of 
the right-hand side of Eq. dnni), we make branch cuts on the a-plane. There are the 
branch points at 01 ^ 2 , the roots of ( 0 , and at ay, and ab±, the roots of C^^\a) = —e 
They are listed in Table 9 as a function of ny. We make branch cuts along the inhnite 
intervals of the real axis {a : —00 < a < ay} and {a : ay, < a < cxo}. There is a 
double branch cut along the interval {a:a 2 <a<cx)}. We also make branch cuts 
along the straight lines joining the points ay,± with the point ay. Now we may continue 
analytically the function (|S3) and determine its boundary values along the sides of the 
right hand cut using the rules of counter-clockwise continuity 0 . Choosing the ways 
along the sides of this cut, we take the limit a ^ 00 . Thus we hnd a pair of complex 
conjugate singular points 


= lim F-\a±iQ) = -hy\l + W^,,{z^))-\ (66) 

a—>00 

where z± = lima^oo C® (® *0), 

z± = {eby)~^ exp(Mi ln(|ai|)/6i - M 2 ln(a 2 )/&i) exp{^mMy/by), 

with My ^2 = (fli — a 2 )~^(l + byay^ 2 )~^- It is obvious from Cauchy’s theorem that if we 
calculate the limits choosing the directions along the left-hand cut, we shall reproduce 
the same values, i.e. 0s,2± = ds.i±- Comparing the modulus of the numbers 9o = — 
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Table 10. The locations of the singularities of a = F{9) and the convergence radii of 
the series (1^ at 3-loops in the theoretical cases when rif = {6 — 16}. 


nf 

IF±„ 

Oo 


P 

6 

IF±i 

1.885 

2.114 

1.885 

7 

IF±i 

3.008 

0.664 

0.664 

8 

IF±i9 

48.17 

0.417 

0.417 

9 

IF^i 

2.08 

0.291 

0.291 

10 

IF^i 

0.761 

0.208 

0.208 

11 

W^i 

0.360 

0.148 

0.148 

12 

W^i 

0.180 

0.102 

0.102 

13 

W^i 

0.087 

0.067 

0.067 

14 

W^i 

0.037 

0.039 

0.037 

15 

W^i 

0.011 

0.018 

0.011 

16 


0.001 

0.003 

0.001 


and 0 s. i± (see Table 10), we determine the radii of convergence of the series (I42|l for 
rif = {6 — 8}. In this Table, we also indicate the reqnired branch of the W fnnction 
which is nsed in formnla ()00|1 . 

Next consider the cases with rif = {9 — 16} (6i < 0, 62 < 0). Then the eqnation 
= —e~^ has a real negative root (oi < Ob < 0) and a pair of complex 
conjngate roots ab± (see Table 9). The solntion to (ITfll which takes real valnes inside 
the real interval at < a < 02 is then given by 

j -b-^{l + Wo{z))-^ if0<a<a2 , . 

—6]"^(1 + W-i{z))~^ if Of, < a < 0, 

where 2 ; = = (e|6i|)“^ exp(T(a)/|6i|). Now we choose the cnts along the real 

axis {a : —00 < a < Oi}, {a : — cx) < a < at} and {a : 02 < a < cx)}. We also 
make cnts along the straight lines joining the branch points ab± with 02. By means 
of the same procednre that was nsed in the previons case we calcnlate the locations of 
the singnlarities = 0s.2±- They are determined by the same formnla (I66|l . The 
relevant branches of W fnnction are listed in Table 10. In this Table, we tabnlated the 
magnitndes of 9q, |0s,i±| = |0s.2±| and p, the convergence radius of the series (jl2]), for 

Uf = {6 — 16}. 


5. The Momentum Scale Associated With the Convergence Radius of the 
Series 

The convergence region of the series (m in the real momentum squared space may be 
easily determined, since the mapping ^ 0 = for real positive > Q\ > 0 is 

monotonic {Q\ being the real Landau singularity of the 2-loop coupling which presents 
if 0 < rif < 8.05). First, we consider the series for large (mainly unphysical) rif values. 
Note that the quantity 0o = f^e Banks-Zaks domain (uf > 8.05) is the infrared 

hxed point of the 2-loop coupling 0 = so that 0 < 0 < for all G (0, cxo). 

From the Table 10, we see that p = 60 inside the interval rif = {14 — 16}. This 
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Table 11. The ratios \/v^^ = and 

the 3-and 4-loop orders for rzf = {0 — 6}. 


rif 

Y ^min 

\/ 

Y ^min 


0 

1.571 

1.525 

1.790 

1.790 

1 

1.566 

1.521 

1.788 

1.788 

2 

1.558 

1.514 

1.784 

1.784 

3 

1.541 

1.500 

1.773 

1.773 

4 

1.505 

1.467 

1.745 

1.745 

5 

1.416 

1.384 

1.678 

1.678 

6 

1.283 

- 

1.623 

1.507 


Qrhp/^ in the MS scheme to 


means that the series (P|) at 3-loops for rif = {14 — 16} converges in the whole interval 
G (0, cx)). Let n}* be the lowest value of n/ for which this equality holds (n}* = 14 in 
the MS scheme). For Uf < the series (j42jl converges in the more restricted domain 
Qmin < < oo (Qmin > 0)' The value of Qmin determined from the equation 

9 = a^‘^\u) = + Wn{zQ))~^ = p ( 68 ) 

where zq = —and u = Q'^/IsA (see Eq. (j2I)). Solving (IHH]) . we obtain 

Mmin = QminM^ = (^1 + exp(p"^). 

The results for the dimensionless quantity = Qmin/A {Qmin = \J Q ^in ) fo the 3- 

and 4-loop orders, for Uf = (0 — 6}, are tabulated in Table 11. We compare with 

the ratio = Qrhp/A, where Qrhp = \J and is the real space-like Landau 

singularity of the coupling. It is seen from the Table that in general the quantity Q^in 
can not be identihed with the real Landau singularity Q^hp- The equality Q^in = Qrhp 
occurs only in the cases where the convergence radius p is determined via the real (space- 
like) Landau singularity. This happens, for example, in the MS scheme in the 4-loop 
case for Uf = (0 — 5}. However, in the case where p is determined via the complex 
Landau singularities, Qi\p±, the relation between Qihp± more complicated. 

Then we observe the inequality Q^in > Qrhp- Such a situation occurs, for instance, in 
the MS scheme in the 3-loop case for Uf = (0 — 5}. 

It is reasonable to compare Qmin with the infrared boundary of QCD, the 
momentum scale pc that separates the perturbative and nonperturbative regimes of 
the theory in the conhning phase. Several estimates for this quantity was suggested 
using different nonperturbative methods. In recent work 1201 an useful nonperturbative 
approximation for the QCD /d-function has been constructed. The model gives a number 
consistent results for various nonperturbative quantities. In particular, it was obtained 
that {pc/A qcd) 3 -Ioop ~ 3.204 and {pc/A qcb) 4 -ioop ~ 3.526 with the perturbative MS 
scheme component of the total /^-function in the 3-th and 4-th orders. Other approach 
to determine the infrared boundary is to use arguments based on dynamical chiral 
symmetry breaking in QCD. There are the results obtained within the nonperturbative 
framework of Schwinger-Dyson equations 1^ . It was found in ref. jl2] that the 
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critical value of the coupling needed to generate the chiral condensate is Oc = vr/4 
(for Nc = 3 QCD). It is reasonable to identify the corresponding scale with the infrared 
boundary 113. To obtain approximations to pc, we may use the perturbative expressions 
for the coupling in the MS scheme. With this simplifying assumption, the equation 
= 7r/4 in the 3- and 4-loop orders, for Uf = 3, yields the estimates 

(/^c/t^Ms)3—loop 1.972 and (/^c/.^ms)4 —loop 2.115. 

The two different estimates considered above are consistent with the inequality 

Thus it seems reasonable to believe that the series expansion dH in the MS scheme 

may be safely used in the whole perturbative region < cx). 

6. Application to Analytic Perturbation Theory 

In the Analytic Perturbation Theory (APT) approach of Shirkov and Solovtsov, 
Euclidean and Minkowskian QCD observables (which depend on the single scale) 
are represented by asymptotic expansions over non-power sets of specihc functions 
and {2t.[f^(s)}))T]^ respectively (see refs. [IH (121), here u = and 

s = s/A^. These sets are constructed via the integral representations in the following 
way 



where the spectral densities to the k-th order are determined from powers of the running 
coupling: = — A(a^^^"'(—-t- *0)). In APT the power series (1^ give rise to the 

following series of functions 


CO 


n=m 

(70) 

00 

n=m 

(71) 

00 

(‘i) = (‘^) m = 1, 2 . . . , 

n=m 

(72) 

where = 1. The sets of coefficients m = 1, 2.. 

., are constructed from 


the set of coefficients of the original series, according to the rules for products 

of power series: \ The spectral densities at the 

2-loop order can be expressed analytically in closed form jsaEni 

with exp (-*7r(l/6i - 1)). (73) 

Now we are going to prove that the series of functions dZOD, (ED) and jza) are uniformly 
convergent over whole ranges of the corresponding variables: 0<m<cxo,0<s<cxo 
and 0 < ^ < cxo. Suppose that the series (17^ is uniformly convergent. Then the series 
dzni) and (ED will also be uniformly convergent. To see this, let us insert the series dZl 
into integral representations given in dnnD and integrate term-by-term. This yields the 
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Table 12. The quantity ^max versus the convergence radii 


nf 

0 

1 

2 

3 

4 

5 

6 

^max 

0.237 

0.237 

0.238 

0.240 

0.243 

0.249 

0.259 

pi3) 

0.594 

0.601 

0.618 

0.653 

0.731 

0.956 

1.885 

p(4) 

0.485 

0.488 

0.497 

0.516 

0.554 

0.641 

0.779 


series (HOD and (HU, which must be uniformly convergent, as results of term-by-term 
integration of the uniformly convergent series. Evidently, the factors l/(<^ -|- m) and l/<^ 
inside the integrals will not destroy this statement. 

Let us now write W_i{z^) = W = X + ly, {I + >V)“^ = 7^exp(*d^), where 
n = {{X + 1)2 + 3;2)-i/ 2 ^ = arcsin(—3^7^) (for the branch hh_i, we have 

—Stt < 37 < 0). According to this, we may rewrite the 2-loop spectral densities as 

= ('^/^i)''sin(n\l'), n = l,2.... (74) 

It is seen from Eq. dZl that the modulus of the spectral densities are bounded above 

l£>f(c)l < (»„»)". (75) 

where ^max = 77- ma x/&i and TT-max is the maximal value of IZ in the range 0 < <^ < cx). We 
hnd useful to use the “Maple 7” gg for determining TT-max numerically. In Table 12, we 
listed numerical values of 6' max in the phenomenologically interesting case Uf = {0 — 6}. 
Note that all the power series ^^. 776 *"', m = 1,2..., have the same radius of 

convergence, as the original series This follows from the dehnition 

( 00 

1=1 

Consider now the set of numerical series of positive terms 

C50 

m = 1.2.... (77) 

n=m 

looking at the numbers in Table 12, we see that ^max is inside the convergence disk of 
the series (|7U|). 0 < ^'max < Hence all the numerical series (|77|) are convergent. 
Combining this fact with the bounding conditions dZU, we hnd that the series of 
functions m = 1, 2 ..., are uniformly convergent by the comparison 

test due to Weierstrass. Then all the series (1Z2D are uniformly convergent. Hence by 
the arguments given above, the series of functions (|7n|) and dZH) are also uniformly 
convergent. 

The series dzni) and ca enable us to calculate the infrared limits of the APT 
expansion functions. Thus we may reproduce remarkable results of Shirkov and 
Solovtsov P Eil in a mathematically rigorous way It is seen from the dehnition 
(p1|l that \im.u^Q+ An\u) = lim 5 ^o+2ti^^(s). Therefore, we shall consider only the 
Minkowskian set of functions. In the 2-loop order, one may calculate the infrared limits 

^ An alternative derivation of these results in the context of the asymptotic solutions to the RG 
equation was recently given in ref. in 


E 
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of the expansion functions using the explicit formulas obtained in ref. [30] • The hrst 
two Minkowskian functions are given by 

= 1 - 7r-i91nW^i(z,), = 7r-i6riSln{iyi(z,)/(l + W,{z,))}, (78) 

where Zg = {ebi)~^s~^^^^ exp{nr{l/bi — 1)). The functions with higher values of index 
are determined by the recurrence relation (see ref. H) 

S'h2(s) = -iV (a?li(5) + ■ (79) 

The asymptotics of Wi{zs) as —>• oo may be determined using Eq. with 
Li = \n.Zs + Ain j^. Combining dZHl) and (|7U|) and taking the limit s —>■ O’*", we hnd 

21(f) (s) ~ + (1 + &i) In-" s + 0(ln-"-i s) ^ (80) 

/Q\ ( 0 \ 

hence lim„^o+( m) = lims^o+2tn (s) = 5^,1. These relations may be extended to 
higher orders by means of the expansions dzni) and m- Thus we can write 


OO 



n=m 


the calculation of the limit of the sum of the series term-by-term, as m —>■ 0, is justihed 
by the uniform convergence of the series. It should be stressed that these results (in 
particular, hniteness of ^^^^(0) originally obtained in [7j) are direct consequences of the 
asymptotic freedom (AF) of the theory. This interesting relationship has been recently 
elucidated by the author of ref. ra using a different technique. In this connection, we 
remark that the recurrence formula dzni) follows directly from the AF, as it was shown 
in refs. isniES- The universality of (0) and 21^^) (0) (the scheme independence and 
invariance with respect to higher-loop corrections) is evident. 

7. Conclusion 

The main objective of this investigation was to study convergence properties of the new 
expansion ©. In section 2 we have systematically discussed the analyticity structure 
of the modihed coupling /IsA) at 3- and 4-loops in the complex u = /IsA plane 

for all rif values in the range 0 < rij- < 16. For higher values of n/, when the (3- 
function has only real zeros, we have reproduced the part of results of refs. HHIIIZI using 
a different technique. For low values of n/, when the /3-function has complex zeros, 
we have determined the analytical continuation of the function t = T{a) choosing the 
cuts properly in the complex coupling plane. With this choice, we have found that the 
running coupling has a pair of complex conjugate singular points in the hrst Riemann 
sheet of the plane besides the real singularity on the positive semi-axis. In many 
cases, just these complex singularities determine the radius of convergence of the series 
© (e.g. 0 < n/ < 5 at 3 loops). 

In section 3 we have proved that in the MS-like schemes the power series 0 has 
a hnite radius of convergence to all orders in perturbation theory for all nj = 0 — 16. 
Therefore, the series inside its circle of convergence represents the exact solution to 
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the RG equation (HH). In the proof, we have used methods of analytical theory of 
differential equations. 

In section 4 we have determined the analytical structure of the modihed coupling 
in higher orders as a function of the 2-loop order coupling 9 {9 = /AA)). We have 

considered the 3- and 4-loop cases for 0 < n/ < 16 and 0 < rij < 7 respectively. We have 
found helpful Eq. (I48j) . the implicit solution for the higher order coupling determined via 
the Lambert-W function. By means of this formula, we have determined the analytical 
continuation of the inverse function 9 = F~^{a) in the complex a plane. This enabled 
us to hnd the locations of the singularities of the coupling a = F{9) in the ^-plane (see 
Tables 6 and 7). The correspondence between the singularities of the coupling in the 
and 9 planes has been established. Comparing various singularities of the coupling in 
the ^-plane, we have determined the radii of convergence of the series (m p(k)_ Pi-om a 
practical viewpoint, the radii of convergence of the original series (jSl), is 

found to be sufficiently large. For example, = 0.965 and p^^'^ = 0.720 at nj = 3. The 
obtained predictions for the convergence radius have been examined by the independent 
numerical calculation. One further important property of the series is high convergence 
rate. In previous papers [S21EB], we observed that in the 3- and 4-loop orders partial 
sums of these series with the first few terms give very good approximations to the 
coupling even in the infrared region. This was confirmed in conventional perturbation 
theory as well as in APT. 

In section 5 we have determined the convergence region of the series (I42|l in the 
momentum squared space. For sufficiently large Uf values (rif > 14 in the MS scheme), 
we have found that the series converges in the whole physical range 0 < < cc. For 

lower Uf values, we have evaluated the lower boundary of the convergence region Q^in- 
We have compared this scale with the estimations of the infrared boundary of QCD, 
obtained within two different non-perturbative approaches and found that < he- 
This is in agreement with the possibility that the series © in the MS scheme may be 
used safely in the whole perturbative region < 0. 

In section 6 we have studied the convergence properties of the non-power series 
constructed from the series (j42|l according to the rules of the QCD Analytic Perturbation 
Theory of Shirkov and Solovstov both in the space- and time-like regions. We have shown 
that the Euclidean and Minkowskian variants of these non-power series are uniformly 
convergent over whole domains of the corresponding momentum squared variables. A 
mathematically rigorous proof of an interesting result of ref. |Zj, the finiteness and 
universality of the analytic coupling at zero momentum, has been presented. 

The series solution 0 may be useful in high-precision calculations of QCD 
observables beyond the 2-loop order in low-momentum regime. It clearly provides 
more accurate results than the standard asymptotic expansion (pH|l for the coupling 
(see ref. [SE!)- This series may be used in different variants of the analytic approach 
to perturbative QCD suggested in refs. jg [m 1221221- It may also be applied in the 
contexts of the “contour improved” perturbation theory of refs. [26-31] and resummation 
methods proposed in refs. El and |22I. Another possible application of the series is to 
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construct the running coupling with consistent matching conditions at quark thresholds 
in MS-like renormalization schemes Ha 113 EDI El. 
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Appendix 


The RG equation to the k-th order reads 

= - E/3n{a=(^^")}”+^ (A.l) 

^ n=0 

the running coupling satishes the normalization condition = ^'^/(dvr), where /i is 

the renormalization point and g is the gauge coupling of QGD. In the class of schemes 
where the /3-function is mass independent /3o and /3i are universal 

/3o = (47r)-^(ll - 2nf/3), pi = (47r)-2(102 - 38nf/3). (A.2) 

The results for the coefficients P 2 and P^ in the MS scheme can be found in refs. jDl] 

and EDI 

P^ = (47r)-3(2857/2 - 5033nf/18 + 325n^/54) (A.3) 


p 3 = (dvr) 
/50065 


_4 /149753 






6 

6472 


162 


+ I n, + 


+ 3564(3 - 
1093 


(1078361 6508, 

+ -T^C3 


V 162 


27 


rif 


729 


-n} 


(A.4) 


here ( is the Riemann zeta-function ((3 = 1.202056903...). The values of the first three 
coefficients 61 ^ 2,3 (bn = Pn/Po^^) in the MS scheme are tabulated in Table A.l. 
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